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Line  - Item  001a 


1 . Causes  of  Ms/m^  Variations  within  a Central-Aslan  Earthquake 
Sequence 

A foreshock-mainshock-aftershock  sequence  in  the  Kirgiz- 
Sinkiang  border  region  is  studied  to  determine  what  factors  other 
than  source  depth  and  focal  mechanism  control  the  surface-wave 
magnitude  of  small  events.  Using  a comparative  event  technique, 
the  mechanism  and  depth  of  the  events  are  deduced  from  azimuthal 
radiation  patterns  of  Love  and  Rayleigh  waves.  All  the  events 
appear  to  be  5-10  km  deep,  with  predominantly  thrust  motion  on  a 
fault  plane  dipping  about  30®  from  horizontal.  Differences  in 
mechanism  and  depth  do  not  adequately  explain  the  dramatic  vari- 
ation in  surface  wave  amplitude  within  the  series  of  aftershocks 
(see  Figure  1).  For  example,  one  aftershock  with  m^  5.2  generated 
surface  waves  roughly  12  times  larger  than  another  aftershock  with 
mj)  5.0. 


The  earthquakes  can  be  divided  into  two  categories:  those 

which  excite  long-period  signals  efficiently,  such  as  the  main- 
shcck  and  several  aftershocks;  and  those  which  do  not,  such  as 
the  foreshock  and  the  majority  of  the  aftershocks.  The  primary 
difference  between  the  two  groups  may  be  the  finite  size  oi*  dura- 
tion of  the  source.  Vfhen  observed  teleseismlcally , the  second 
group  typically  generates  a very  simple,  short-period  P-wave  con- 
centrated within  the  first  1 to  2 s (upper  trace.  Figure  2),  sug- 
gesting nearly  a point  source.  However,  earthquakes  within  the 
first  group  generate  a much  more  complex  short-period  P-wave, 
with  substantial  energy  spread  out  over  the  first  10  to  15  s,  as 
shown  in  the  lower  trace  of  figure  2 (duration  of  time  marks  is 
2 s).  The  maximum  amplitude  of  the  short-period  signal  is  not  a 
very  complete  description  of  its  character.  A better  description 
would  include  both  amplitude  and  duration,  as  when  the  area  of 
the  envelope  of  the  short-period  record  is  measured.  A comparison 
of  these  two  measures  of  earthquake  size  is  given  in  Figures  3 and 
4,  which  show  log  of  maximum  amplitude  (i.e.,  mp ) and  area  of  short 
period  envelope,  respectively,  when  plotted  versus  Mg.  Using  the 
area  of  the  envelope  in  place  of  mp  will  remove  much  of  the  scatter 
in  Ms/mp  relations,  without  greatly  increasing  the  effort  involved 
in  making  the  measurements.  These  observations  suggest  that  finite 
ness  of  the  source  can  play  an  important  role  in  controlling  Ms/mp, 
even  for  earthquakes  with  body  wave  magnitude  4.8  or  smaller. 

2 . Distribution  of  Great  Shallow  Earthquakes 

Based  on  a world-wide  survey  of  great  earthquakes,  Kelleher 
and  McCann  (1976)  found  that  the  distribution  of  large  shallow 
earthquakes  along  subduction  boundaries  does  not  agree  with  the 
distribution  pattern  that  might  be  predicted  from  a simple  model 
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of  plate  tectonics.  That  Is,  alonp  extensive  sections  of  some 
island  arcs  large  shocks  occurred  Infrequently  or  not  at  all 
during  recorded  history.  Most  of  these  zones  of  long-term  qui- 
escence are  nearly  coterminous  with  segments  of  the  margin  where 
zones  of  seamounts,  aseismic  ridges  or  other  bathymetric  highs 
of  the  underthrust  slab  appear  to  be  interacting  with  the  subduc- 
tion  process.  This  spatial  correlation  suggests  that  at  least 
some  of  the  long-term  absences  of  great  shocks  may  result  from 
a tectonic  origin  and  not  from  temporary  intervals  of  strain 
accumulation.  Thus,  major  departures  from  classic  subduction 
activity  may  develop  where  significant  bathymetric  features  inter- 
act with  a convergent  margin. 

Prom  a similar  survey  of  the  source  regions  of  seismic  sea 
waves,  or  tsunamis,  Kelleher  (1976)  Indicates  that  many  have 
Important  features  in  common,  and  that  at  the  present  time,  there 
exists  a number  of  sites  which  appear  to  fit  these  occurrence 
patterns  of  past  tsunamis.  Probably  a majority  of  great  tsunamis 
were  generated  along  subducting  margins  by  low-angle,  thrust-type 
earthquakes  which  had  rupture  zones  extending  as  much  as  several 
hundred  kilometers  along  strike.  For  these  great  earthquakes, 
therefore,  large  vertical  deformations  of  the  sea  floor  occurred 
because  of  the  large  seismic  moments  (proportional  to  fault  area 
X displacement),  and  despite  the  typically  low  angles  of  faulting 
(10°  to  30°).  Crude  estimates  of  the  size  and  location  of  some 
future  shocks  of  this  important  class  of  earthquakes  are  possible 
by  the  technique  of  seismic  gaps  (segments  of  major  plate  boun- 
daries that  have  not  ruptured  for  many  years).  Thus,  significant 
seismic  gaps  along  subduction  boundaries  deserve  attention  as 
potential  locations  for  future  tsunamis. 

A major  new  emphasis  in  research  concerns  the  development  of 
the  inner  wall  of  the  trench  and  its  relationship  to  the  size 
and  frequency  of  great  earthquakes.  The  presence  of  numerous 
scarps  and  large  deep-sea  terraces  within  a well-developed  pattern 
of  imbricate  thrust-faulting  may  be  intimately  related  to  and 
possibly  an  indicator  of  the  episodic  release  of  tectonic  strain 
through  great  earthquakes. 


3.  Volcanic  Structures  as  Indicators  of  Tectonic  Stress 
Orientation 

A new  method  for  obtaining  tectonic  stress  orientation  from 
volcanic  structures,  proposed  by  Nakamura  (1969  and  in  press), 
was  applied  to  the  Aleutian  and  Alaskan  volcanoes  and  volcanic 
fields.  The  method  is  essentially  the  recognition  of  the  pre- 


ferred orientation  of  radial  and  parallel  dike  swarm  development 
by  means  of  their  probable  surface  manifestations,  such  as  flank 
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volcanoes  and  volcanic  fleldn  primarily  the  trend  of  maximum 

compression  of  the  horizontal  components  of  the  tectonic  stress. 
When  the  method  is  applied  on  a regional  scale,  one  can  further 
identify  the  trend,  either  as  the  maximum  or  as  the  intermediate 
axis  of  the  tectonic  stress. 

/ 

The  trends  of  the  volcanic  features  which  Indicate  the  maxi- 
mum compression  of. the  horizontal  components  were  obtained  from 
fifteen  volcanoes,  including  Buldlr,  the  westernmost  Aleutian 
volcano,  and  Illamna  volcano,  near  the  eastern  end  of  the  Aleutian 
volcanic  belt.  Trends  of  the  volcanic  features  generally  coincide 
well  with  the  azimuths  of  slip  vectors  for  the  relative  motion 
between  the  Pacific  and  North  American  plates  (Minster  e^  , 
197^),  thus  providing  evidence  that  the  obtained  horizontal  direc- 
tions at  these  volcanoes  are  for  the  maximum  compresslonal  axis 
of  the  tectonic  stress. 

General  east-west  trends  were  obtained  from  nine  volcanoes 
and  volcanic  fields  on  Islands  and  at  the  mainland  coast  of  the 
Alaska-Bering  Sea  shelf.  These  volcanoes  and  volcanic  fields 
are  mostly  of  alkali-basalt  (Smith  et  al . , 1973)  and  are  generally 
associated  with  normal  faults  of  slinTlar  trends.  Therefore,  the 
obtained  directions  are  most  probably  those  of  the  intermediate 
axis,  with  the  maximum  compression  axis  in  the  vertical  direction. 

These  results  give  strong  support  to  the  idea  that  the  flank 
eruption  on  polygenetic  volcanoes  can  be  regarded  as  large-scale 
natural  experiments  for  magmatic  hydrofracturing.  Moreover,  they 
have  important  implications  for  the  tectonics  of  island  arcs  and 
back-arc  regions  (such  as  the  Kuril-Kamchatka  arc  and  the  Sea  of 
Okhotsk  in  the  U.S.S.R.):  (1)  Volcanic  belts  of  some  Island  arcs, 

at  least  the  Aleutian  arc,  are  under  compresslonal  stress.  For 
such  arcs,  it  is  Improbable  to  separate  the  various  parts  of  the 
arc  and  form  a marginal  basin;  (2)  The  compresslonal  stress  at 
the  arc,  probably  generated  by  the  underthrusting,  appears  to  be 
transmitted  across  the  entire  arc  structure,  but  is  apparently 
replaced  within  several  hundred  kilometers  by  a different,  ten- 
slonal  stress  system  in  the  back-arc  region;  (3)  Both  stress 
systems,  the  compresslonal  and  the  tensional  ones,  are  associated 
with  volcanoes  and  volcanic  fields,  which,  however,  differ  in  the 
chemistry  of  their  magma.  This  difference  supports  the  idea  that 
the  back-arc  stress  system  has  its  own  source  at  considerable 
depth  beneath  the  crust. 

^ . Volcanic  Eruptions  and  Their  Role  in  Seismic  Gap  Theory 

Volcanic  eruptions,  now  perhaps  predictable  events,  may  be 
used  in  long-term  earthquake  prediction.  Klmura  and  McCann 
(in  preparation)  examine  the  possible  relationship  between  major 


volcanic  eruptions  and  the  occurrence  of  shallow  thrust  earthquakes 
alonp,  the  I’aclflc  marf’iln.  Althour.h  this  study  Is  still  In  Its  nre- 
llmlnary  sta^on,  none  noteworthy  results  can  now  he  reported. 

Volcanic  activity  alonp  the  Kamchatka-central  Kuril  suoductlon 
zone  was  relatively  active  during  the  period  1923  - 19^0.  During 
the  period  19^0  - 1952  only  five  volcanoes  v/ere  active,  three  in 
north  Kamchatka  and  two  In  the  Kuriles.  In  1952  the  great  Kamchatka 
earthquake  (Mg  8.25)  occurred,  rupturing  500  km  along  southern 
Kamchatka  and  the  northern  Kuriles.  After  this  event  volcanic 
activity  increased  north  and  south  of  the  rupture  area  but  remained 
at  a low  level  in  the  region  near  the  rupture  zone.  Similar  rela- 
tionships between  volcanic  eruptions  and  the  large  thrust  earth- 
quakes can  be  found  (e.g.,  1952  Aleutians,  I960  Chile,  1968  Hok- 
kaido) . 

One  possible  model  to  explain  the  observed  relationship  is  one 
which  considers  the  elastic  nature  of  the  plate  boundary  in  re- 
sponse to  stresses  applied  to  it.  Prior  to  a large  thrust  event 
stresses  would  increase  on  a plate  boundary.  This  stress  would 
be  regional  in  extent  and  involve  a region  several  times  the  area 
of  any  future  shock.  This  increased  stress  would  induce  volcanic 
eruptions  through  compression  of  magma  source  regions.  If  the 
future  event  is  regional  in  extent  (e.g.,  rupture  zones  several 
100' s of  km  large)  then  the  decade  of  decreasing  activity  before 
the  events  could  be  due  to  dilatancy  phenomena  reducing  the  pore 
pressure  in  the  magma  source  region  and  therefore  the  need  for  a 
volcanic  eruption.  After  the  earthquake  stress  in  adjoining  re- 
gions (gaps)  would  markedly  increase  especially  in  those  regions 
that  had  been  free  of  large  thrust  events  in  the  several  previous 
decades.  This  increase  in  tectonic  stress  would  be  reflected  in 
an  increase  in  volcanic  activity. 

Line  - Item  001b 


1 . A Low  Velocity  Layer  over  the  Tibetan  Plateau 

Surface  wave  data  from  more  than  ten  Tibetan  earthquakes  have 
been  digitized  and  moving  window  analyzed.  The  period  range  is 
5 - 100  seconds  with  the  corresponding  group  velocity  window  of 
2.3  - ^4.1  km/sec.  The  observed  dispersion  curves  have  been  made 
and  compared  with  their  respective  theoretical  curves. 

We  have  constructed  for  the  Tibetan  Plateau,  a number  of  crustal 
and  upper  mantle  models  which  we  consider  represent  excellent  aver- 
age structures  along  the  travel  paths.  One  of  them,  TP-3  model, 
is  shown  in  the  following  table; 


RAYLEIOH  WAVE 


MODEL  TP- 3 


Layer 

No. 

Thickness 

Dep . 

P 

Velocities 

S 

Velocities 

Density 

1 

3.5000 

3.5000 

4.5000 

2.6000 

2.4000 

2 

34.5000 

38.0000 

5.9800 

3.4500 

2.8000 

3 

30.0000 

68.0000 

6.3000 

3.6400 

2.9000 

4 

INF 

INF 

8.0000 

4.6000 

3.4000 

The  main  findings  contained  in  our  preliminary  report  still 
hold  true  except  that  there  now  appears  to  be  some  doubts  regard- 
ing the  existence  of  a relatively  undeformed  layer  of  low  velocity 
sediments.  Since  the  existence  (or  non-existence)  of  the  unde- 
formed layer  of  low  velocity  sediments  proposed  by  Chen  and  Molnar 
(1975)  has  Important  tectonic  implications,  we  have  made  some  pre- 
liminary investigation  into  it.  We  think  that  the  surface  wave 
dispersion  across  the  Ganges (to  reach  HDI  and  LAH  stations)  and 
the  complexities  in  the  Himalayas  have  such  significant  effect 
(especially  in  the  short  period  range)  that  care  has  to  be  fully 
exercised  when  analyzing  the  high  frequency  surface  waves. 

Line  - Item  001c 

1.  A Study  of  the  Detection  Capabilities  of  Several  of  HGLP 
Seismograph  Stations 

For  the  first  study  a swarm  of  earthquakes  that  occurred  on 
the  mid-Atlantic  ridge  in  May  of  197^  near  27. ^®N  and  44.4®W  is 
being  used  as  a test  group  to  determine  how  accurately  a few 
HGLP  stations  can  locate  an  earthquake  and  determine  its  focal 
mechanism.  This  swarm  was  selected  because  it  is  within  50°  of 
three  HGLP  stations:  Ogdensburg,  N.J.,  Kongsberg,  Norway,  and 

Toledo,  Spain.  The  location  (the  distance  and  azimuth)  of  each 
earthquake  in  the  swarm  will  be  recomputed  using  the  three  HGLP 
stations  and  compared  to  the  PDE  locations  as  the  standard.  In 
addition  to  locating  the  events,  magnitudes  and  focal  mechanisms 
will  be  derived  using  data  from  only  the  three  HGLP  stations. 

These  will  be  compared  with  those  obtained  by  normal  multi- 
station methods.  Several  digital  filtering  techniques  have  been 
employed  to  Improve  the  detection  capabilities  of  the  HGLP  sta- 
tions. Band-pass  and  master  event  are  among  the  techniques  being 
used. 
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2 . Possible  Fine  Structure  In  Underthrust  j np,  Lithosphere 

Many  studies  have  been  reported  concernlnfr  the  development  of 
trench-arc-backarc  systems.  Some  of  them,  however,  are  quite 
Ignorant  of  recent  geophysical  observations  In  these  very  complex 
regions.  Dr.  Yoshll  has  compiled  some  geophysical  observations 
In  the  northeastern  Japan,  one  of  the  most  tynlcal  trench-arc- 
backarc  systems.  From  the  abundant  observation  data,  maps  and 
cross  sections  for  bathymetry,  crustal  structure,  gravity  anomaly, 
heat  flow  and  earthquakes  were  constructed. 

The  extensive  analyses  of  earthquakes,  gravity  and  heat  flow 
have  been  made  so  far.  Reconstruction  of  the  earthquake  hypo- 
central  cross  section  by  using  pP-depths  reported  by  the  Inter- 
national Selsmologlcal  Centre  Implies  surprising  concentration 
of  the  earthquakes  on  a nlane  of  only  15  km  thick  extending  from 
the  trench  axis  down  to  the  depth  of  about  300  km.  Moreover,  the 
second  plane  can  be  clearly  seen  about  ^0  km  below  this  main 
plane.  Focal  mechanism  solutions  of  several  events  on  this  highly 
detailed  earthquake  distribution  help  our  understanding  of  complex 
motions  acting  on  the  trench-arc-backarc  system,  such  as  bending 
and  unbending  of  the  oceanic  plate  and  thrustal  motion  between 
the  continental  and  oceanic  plates. 

Line  - Item  OOlf 

1 . Numerical  Techniques  for  Computing  Seismograms  of  Several 
Body  Phases 

A study  Is  In  progress  of  a technique  that  calculates  seis- 
mograms for  body  waves  which  undergo  turning  point  Interactions. 
The  method  Involves  numerical  (but  cheap)  computations  In  the 
complex  ray  parameter  plane,  of  a type  which  has  Important  dif- 
ferences from  Cagniard's  method.  We  use  a complex  velocity  pro- 
file to  account  for  attenuation,  and  the  uniformly  asymptotic 
method  of  Langer  to  account  for  frequency-dependence  of  reflec- 
tion-transmission coefficients.  The  ray  systematics  for  studying 
a high  velocity  region  over  a low  velocity  (e.g.,  the  core-mantle 
boundary),  with  associated  tunnelling  effects,  have  been  known 
for  several  years.  However,  the  (geophysically)  more  common  case, 
of  low  velocity  region  over  high,  leads  to  ray  systematics  which 
have  only  recently  been  understood.  The  work  of  Nussenzvelg 
which  correctly  evaluates  all  the  waves  In  the  high  velocity 
medium  which  are  repeatedly  refracted  back  up  to,  and  reflected 
from,  the  low  velv^city  medium.  The  Debye  ray  expansion  Is  not 
made.  In  the  region  where  it  is  poorly  convergent. 

Initial  Indications,  from  seismograms  observed  and  calculated 
near  the  point  D of  PKP  waves  (which  is  the  critical  distance  for 
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K waves  incident  from  fluid  core  upon  solid  inner  core),  are  that 
impedance  Jumps  are  greater  than  had  been  supposed,  and  that  the 
amplitude  maximum  associated  with  the  critical  distance  is  displaced 
to  shorter  distances  at  lower  frequencies.  We  have  also  found  that 
computed  seismograms  for  SmKS  (m  = 1,  2,  ...)  exhibit  many  non- 
ray effects  in  the  distance  range  100°  - 125°.  Such  waves  are 
an  excellent  example  of  the  way  in  which  successively  hlgheir- 
order  reflections  (SKKS,  SKKKS,  etc.)  have  arrivals  within  the 
waveform  of  previous  reflections.  The  total  seismogram  for  all 
the  multiples  is  exceedingly  sensitive  to  velocity  gradient  in  the 
outer  fluid  core,  and  initial  indications  are  that  the  Jeffreys- 
Bullen  gradient  is  correct. 

We  are  currently  testing  our  method,  for  S waves  interacting 
with  upper-mantle  discontinuities. 
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